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NOMENCLATURE 


Coefficients of the velocity distribution 
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(Appendix D) 
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Temperature 
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Velocity component along the body surface 
(Appendices C and D) 
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Bulk or mass averaged velocity 
Nondimensional mass averaged velocity. 
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Nondimensional thermal velocity. 
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Molecular (or atomic) weight of species i 
Molecular weight of mixture 
Coordinate parallel to body (Fig. 1) 

Coordinate perpendicular to body (Fig. 1). 
Coordinate parallel to body (Fig. 1) 

Energy accommodation coefficient 
Body angle defined in Fig. D1 
Kronecker delta 

1/2 

Reynolds number parameter, [u ref /po.Moo r N] 

Recombination coefficient (or, fraction of 
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Ratio of specific heats 

Surface curvature 

Mean free path 
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Perturbation part of distribution function 

General property of ith species such as mass, 
momentum and energy 
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Circumferential angle (Appendices D and E) 
Density 

Accommodation coefficient (or fraction of 
incident particles that stick) 
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Superscripts 


Atom 

Equilibrium value 
Molecule 
Specie indices 
Coordinate indices 
Specie index 
Reference quantity 
Component directions 
Normal component 
Tangential components 
Freestream value 


+ 
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Dimensionless quantity; also used locally 

Incident flux 

Specularly reflected flux 



SURFACE-SLIP EQUATIONS FOR MULTICOMPONENT, 

NONEQUILIBRIUM AIR FLOW* 

Roop N. Gupta** 

Carl 0. Scottt 
and 

James N. Mosstt 
SUMMARY 

Equations are presented for the surface slip (or jump) values of spe- 
cies concentration, pressure, velocity, and temperature in the low-Reynolds- 
number, high-altitude flight regime of a space vehicle. These are obtained 
from closed-form solutions of the mass, momentum, and energy flux equations 
using the Chapman -Enskog velocity distribution function. This function 
represents a solution of the Boltzmann equation in the Navier-Stokes ap- 
proximation. The analysis, obtained for nonequilibrium multicomponent air 
flow, includes the finite-rate surface catalytic recombination and changes 
in the internal energy during reflection from the surface. Expressions for 
the various slip quantities have been obtained in a form which can readily 
be employed in flow-field computations. A consistent set of equations is 
provided for multicomponent, binary, and single species mixtures. Expres- 
sion is also provided for the finite-rate species-concentration boundary 
condition for a multicomponent mixture in absence of slip. 

INTRODUCTION 

For an accurate prediction of the aerothermal environment of a space 
vehicle entering the earth's atmosphere in the high-altitude low-Reynolds- 
number flight regime (Ref. 1), the multicomponent, nonequilibrium gas 
chemistry as well as the wall slip and catalysis effects, must be evaluated. 

*The authors are grateful to Dr. Fred G. Blottner (of Sandia National 
Laboratories, Albuquerque, NM) for many useful discussions pertaining to 
the mass flux expressions contained in Appendix B. 

**01 d Dominion University, Norfolk, VA 
t Johnson Space Center, Houston, Texas 
ttLangley Research Center, Hampton, VA. 



Such flow fields are of particular Interest for aeroassist and space shuttle 
vehicles. This study was undertaken to obtain the boundary relations that 
incorporate the effects of slip, multicomponent diffusion, wall catalycity, 
and changes in internal energy of the molecules (during reflection from the 
surface) for application to flow-field calculations under the general as- 
sumption of local thermodynamic equilibrium. 

Scott (Ref. 2) first presented the wall boundary conditions for a 
multicomponent mixture with diffusion and wall-catalysed atom recombination. 
In obtaining these boundary conditions, he used a first order velocity dis- 
tribution function at the edge of the Knudsen layer next to the wall, where 
the continuum model of the gas breaks down. These boundary conditions, 
obtained from the kinetic theory considerations, provide solutions at the 
top of the Knudsen layer that would match the solution of the Navier-Stokes 
equations in the bulk outer flow. Hendricks (Ref. 3), using Scott's formu- 
lation, obtained simplified expressions for engineering applications with 
some corrections to Scott's expressions. Hendricks' analysis, however, 
contained some gross errors in obtaining the engineering expressions. This 
paper reanalyses the wall boundary equations by using the approach of refer- 
ence 2 and provides appropriate relations for the various quantities with 
surface slip in a form which can readily be employed for flow-field computa- 
tions. An effort has also been made to reconcile the differences between 
slip expressions employed by the different researchers. The present analy- 
sis provides a consistent formulation for the slip equations for a multi- 
component, binary, and single species mixture. Expression is also provided 
for the finite-rate species-concentration boundary condition for a multi- 
component mixture (in addition to that for a binary mixture) in the absence 
of slip. This may be of interest for the shuttle flow-field calculations 
based on multicomponent diffusion (Ref. 4). 
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The main difference (among other details) between the results contained 
in reference 5 and the present work is that the various internal degrees of 
freedom for a molecule were considered frozen during reflection from the 
surface in reference 5, whereas, they are allowed to change in the present 
analysis.*** Further, the results are provided in both spherical and body- 
oriented coordinate system. 

ANALYSIS 

The slip conditions are taken to exist across the Knudsen layer, which 
is on the order of one mean free path in thickness as sketched in Figure 1. 
The analysis outlined here follows the approach of References 2 and 5. It 
is based on Shidlovskiy 1 s (Ref. 6) assumption that the distribution function 
near the wall can be described to first-order accuracy by the so-called 
Navier-Stokes approximation. However, a deviation is made from the proce- 
dures of Shidlovskiy in that a Chapman-Enskog type distribution function for 
a multicomponent mixture obtained by the variational method of Hirschfelder, 
Curtiss, and Bird (Ref. 7) is used. The Chapman-Enskog distribution 
function allows accounting for diffusion. The analysis contains the follow- 
ing assumptions: 

(i) The energy and momentum accommodation coefficients (i.e. a 
and 0, respectively) have the same value. 

(ii) The fluxes of mass, momentum, and energy across the Knudsen 
layer are assumed constant. This is consistent with the 
assumption of negligible variation of the velocity distri- 
bution function through the Knudsen layer. 

***Since the assumption of local thermodynamic equilibrium is employed, only 
those internal energies are considered which equilibrate readily with the 
translational energy. 



(iii) The internal energy associated with the rotational and 
vibrational modes readily equil ibrates with translational 
energy. 

The interaction model at the gas-solid interface, with the various 
fluxes sketched in Figure 1, can be mathematically stated for dissociated 
air as provided through the following equations. 

For a recombining atom : 

F A - F^ + (l-e A ) F a +.(e A -v A ) Fjj ; A * 0,N (1) 


For a molecule gaining from the corresponding atom recombination : 


+ < l -V F n - 6 m f h + 'A r A • M * 


/ ~ V 

u; 


For all other atoms and molecules (surface is assumed to be noncata- 
lytic with respect to them): 

F, • F 1 , + (1-8,) F, +8 1 F i w (3) 

where F.. denotes a convective property such as mass, momentum, or energy. 

Summing over all the species, we obtain from Eq. (1), (2), and (3), the 
following expression for the net flux of momentum or energy: 



That is, the net flux at the outer edge of the Knudsen layer equals the 
incident flux, plus the specularly reflected flux (incident minus the frac- 
tion that sticks) at the wall, plus the diffusely reflected flux (those that 
accommodate to the wall) from the wall. 

Each species is treated separately in the mass balance equations. 
Therefore, Eqs. (1), (2), or (3) are employed depending on the species being 
considered. In Eq. (1), the diffusely reflected flux consists of those 
atoms that are accommodated to the wall minus those that recombine. For the 
molecules in Eq. (2), the diffusely reflected term is present along with the 
source term resulting from the appropriate atoms recombining on the surface. 

VARIOUS FLUXES AND THE DISTRIBUTION FUNCTIONS 
The interaction model of Eqs. (1) through (4) is employed to obtain the 
slip boundary equations at the gas/solid interface. Through Eq.(l), the net 
fluxes of species, momentum, and energy at the outer edge of the Knudsen 
layer are equated to the difference between the incident and reflected 
fluxes at the wall. These fluxes are assumed to be constant across the 
Knudsen layer and are obtained from moments of the distribution function. 

For a convected property <|>^(v) such as mass, momentum, or energy for the 
ith species, the net flux of that property normal to the wall at the outer 
edge of the Knudsen layer, for example, is 

f,-l‘ C C (?) (5) 

where V^' is the normal component of the molecular velocity and f^ is the 
velocity distribution function at the edge of the Knudsen layer. 

Similar integrals are obtained for the incident and reflected fluxes by 
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integrations over appropriate half-spaces in molecular velocity: 
Incident flux: 


F* -£ L C v’ + i (V) fj(?) dBV 1 (6) 

Specularly reflected flux: 

^ * C r 0 C V y ♦ 1 (^)-F s< V x- - V y. V z )d3V 1 (7) 

Oi ff usel y ref 1 ected f 1 ux : 

*?■£/; CvJ* 1 #) '!(*)■ a** 1 (8) 

where is the Maxwellian velocity distribution function. 

n 

The velocity distribution functions used in integrals contained in 
relations (5) through (8) are those for a nonuni form multicomponent mixture 
perturbed out of equilibrium: 


f 1 ^) * f (0)l (1/) [1 + fr 1 (V)] (9) 

where f^^(V) is the Maxwellian distribution function for the ith species 
given as 


f (0)i (tf) 


3/2 


n i m i 

(2trkT) 


3/2 


-(m i /2kT)V. 2 
e ! 


( 10 ) 



and 


«’(to = 


i 3(1 nT ) 

rt — 

,x k 


B^*nf C^J) d J 

■9x t j-i K K 


Here k and l are the dummy indices for three coordinate directions and 
the summations with index j represents summation over all the species. 
The summation convention for repeated indices is used. The coefficients 
A k* C and are functions of the dimensionless velocity: 


and are defined as 



*k = ho + a ii(f- w ?)l w k 


b 10 t w k < -T W 1 


C 1(J) 3 

L k 


c (j) yi 

c i0 w k 


where a^Q, a.p b^ Q , and c^g ' are constants determined from the varia- 
tional problem in the first approximation for a mixture as given in ref. 7 
and is the Kronecker delta such that 


3 if k 3 * 
* kl = 0, if k * i. 
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These constants are functions of the collision integrals and are related to 
the transport properties. The simplified form of these constants is pro- 
vided in Appendix A of this paper. More details can be found in Appendix A 
of Ref. 2, or ip Ref. 7. The diffusion vector d^ in Eq. (11) is related 
to the diffusion velocity of the jth species and is defined after neglecting 
the external forces as (Ref. 7): 


a r n j n 


3x k n 


U) * ( J- 

n 


n.m. . 

L_ (inp) 
p 3xi, 


(16) 


where n. and m. are, respectively, the number density and mass of the 

J J 

jth species, n is the total number density n = T n«, p is the total 

0 J 

mass density p * l n.m., and p is the total pressure p = l p. . A 

j Y j j 

simplified form for d^ is provided in Appendix A. 

The total mass averaged velocity Vq^ (i.e. the kth component) 

appearing in Eq. (11) is defined as 


Ok 


p 


NS 
l n.m 


j=l 


n V k 


where v^ is the mean of total velocity vjj! = Vq^ + of the jth 
species averaged over the distribution function and is the thermal (or 
peculiar) velocity, also introduced in Eq. (5). The thermal velocity Vj| 
of the jth species averaged over the distribution function is known as the 
diffusion velocity V£: 





"j 


On f J rt> *** 


( 17 ) 


A simplified expression for the diffusion velocity in terms of transport 
properties is provided in Appendix A. 

THE BALANCE AND SLIP EQUATIONS 

Based on Eqs. (1) through (4), the balance equations for the ith 
species for fluxes normal to the surface of the species mass M.^, the normal 
component of momentum P^, the tangential component of momentum and 

the energy E^ are obtained as follows: 

(a) Species mass flux 

For a recombining atom 


V ■ e A < + < 9 A - V 1$ S A = 0, N 
For a molecule gaining from the corresponding atom recombination 


(18) 


M My = 9 M i + e M Kl + 7 A ; N . 0 2 , Na 


(19) 


For all other atoms and molecules 


M, a 8 . + 9 . m'? 

iy i l i i 


( 20 ) 


(b) Normal momentum flux 


NS NS , NS 

l p. = y (2-0 .) p. + y 0 . P^ = P 

i=i ^ i*i ' 1 i=i 1 1 y 


( 21 ) 
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(c) Tangential m omen tun flux 

NS NS 

h Pf M ‘, z .i 9iPl n 

t 

(d) Energy flux 


NS 


NS 


E y ' 


y e. * y e, + y e. m 
t-, iy A, iy h iy 


i*l 

NS 


i*l 


diatomic 

molecules 


iy iy 


ns Tl ns T 

l 9 E* ♦ I 8 E, TW + I 9 eft 

i=l 1 1 i=l 11 i 1 1 1 


diatomic 

molecules 

+ 1 w„w 

1 9 i e i M i 
diatomic 
molecules 


( 22 ) 


(23) 


where e^ is the internal energy of ith species that readily equilibrates 

T 

with the translational energy under the assumption of local thermo- 
dynamic equilibrium. For example, 

rotation kT 


In writing the energy flux balance of Eq. (23) it has been assumed that 
there is no change in internal energy during specular reflection. 

In obtaining Eqs. (18) through (23) we have used the following 
relations 

tThe energy balance is based on the assumption that the various energies 
considered readily equilibrate with the translational energy. 



M 


+ 

i 






P'Jll = 0 , and 



Because it is assumed that the atoms are consumed at the wall by catalytic 
recombination in Eq. (18), the net mass flux M. * 0. Similarly, the net 
mass flux M My * 0 in Eq. (19). However, = 0 in Eq. (20) for the 
atoms and molecules for whom the surface is assumed noncatalytic. 

By substituting Eqs. (5), (6), (7) and (8) with the definition of 
( V) as mass for Eqs. (18) through (20), as normal component of momentum 
for Eq. (21), as tangential component of momentun for Eq. (22) and as energy 
for Eq. (23), respectively, and carrying out the integrations*, one obtains 
equations relating the slip properties to wall properties and gradients at 
the edge of Knudsen layer. All accommodation coefficients 0. are assumed 
to be equal to 0 . 

Number density (or concentration) slip (obtained from mass flux balance): 



Pressure slip (obtained from the flux balance of normal component of momen- 
tum) 


tThe mass, momentum, and energy fluxes in terms of evaluated integrals over 
the distribution function are given in Appendix B. 
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„ 3 V_ 3V rt 3 V rt 

p 1- - [■[_*_♦ _£- 2 -*>]. 


3 X 


32 


3y 


NS 

+ [k- I ni(a j 0 -i.a n )] s+ I p w 

m 3y i=l 1 10 2 11 s J 


/{ 1 + £± I [n. f c^ ! d J ] } 

2 vr 1 1 J=1 10 


(25) 


NS 


« I w 

with viscosity u * L kT Y n, b 


2 1-1 


iO 


Velocity slip (obtained from the flux balance of tangential component of 
momentum) 


v 0x 58 1 ^ (— ) »c ( 


3 V- 3v 


Ox 


+ _°2) 


q V / C 

20 * 3y 3x . 


NS 


.S r 3A nT 


+ - I Pit 
2 i=l 1 3x 


( a iO “ ~ a il^ 


NS 


n I ctf dj $} 


J-l 


NS 

ill 


V/% * 

Oz 



(26) 


9v n 

+ _J^ 

3 Z 


. NS _ - NS ... . 

i I P [ii2l(« 10 -ia ) - R I c ) dJ] s } 

? i=l 1 3 z 1U 2 11 j=l 1U z s 
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Temperature slip: 



Equations (24) through (28) differ from the corresponding expressions 
provided in Ref. 2 due to small errors and also due to the differences in 
the interaction model anployed at the gas-solid interface. 

The constants a iQ , a^, b iQ , and c^ (also known as the Sonine ex- 
pansion coefficients) appearing in Eqs. (24) through (28). may now be ex- 
pressed in terms of the transport properties as given in Appendix A. Using 
these relations along with the various flux expressions of Appendix B and 
also expressing d^ in terms of the gradient of mass fractions 8C./3y as 
given in Appendix A (by neglecting diffusion due to pressure gradients and 
external forces), the following equations are obtained after some algebraic 
simplifications:^ 


§The approximations made in the expressions for a^ and b^Q are given in 
Appendix A. 



Number density (or concentration) slip 


Same as Eq. (24). 


Pressure slip 


p* = 1- •_[, ( !I2 l ♦ - 2 !>)] s ♦ i (i±) [I ii f m. 

3 3 x 3 z 3y s 5 iT » 3y i=l 1 1 J 2kT 


*- P W 1/(1*2{^-) I 
2 2 Vi T i = 1 


2-q> v S ri 


DO . 

*J, VI— 


NS S 

Cj 

3 q=i W„ 


Velocity Slip 


3v„. 3 v. 


'S-'Ws 


I 

& "1^ K ll 


I1J IW 

-n s .1 i"i l D ij[— 
1=1 j=l J 3 X 


NS = 3C NS 

c ji ( - -5)] s i n 

J q=l m q 3x s i=l .. 1 y 1 


uz \ 2 9 ^kT ay 


2E. + Qy.)i + 1 r_L !I 

y 3Z s 5 = kT 3 Z 
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Temperature slip 


w 


NS M. 


m.e. w M. 


{ _V¥ T — l ( 1 1 ) iy -V^ (2-8), 

i ' 

's diatomic 
molecules 


- 1=1 m. 

"s 1 


2 n. 


kT w m i 2n s 9 


w s, 


diatomic 

molecules 


w i 


m.e. ^ M. , / 2k T m e (e.-e. )m. 

I (- 7^3 — + 7 l J— (-) C 

A m- m. fcT 


diatomic 

molecules 


i l 


w 



NS M, 


/(-vf (i±> a i ii- * i Ji.) s 

Q 2 p By 4 i=l n m. 


, NS / 2 kT_ m_ P, 

I I / 1 (— ) C, (3 -1* + 1)} 

4 i=l \l m, m. n s 

i i p ^ 


(32) 
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where the mass fraction C i and the mass of a mixture molecule m are 
defined, respectively, as 


and 

NS _! 

m - ( I C./m.) 
j=l J J 

Equations (24), (29), (30), (31), and (32) differ from Eqs. (25) through 
(28) of Reference 3. For dissociated air, the gas-surface interaction model 
employed in Reference 3 appears to be inconsistent. Further, the diffusion 
vector dj^ appears to have been incorrectly evaluated in Reference 3. If 
one carries out the simplifications in Eqs. (23), (24), and (28) of Refer- 
ence 3 (which contain the binary assumption) through the evaluation of 

— (n./n) in terms of mass fractions as given by Eqs. (35), (36), and (39) 
ay J 

of the same reference, erroneous results are obviously obtained because 
NS 

l (ac./ar) (which is zero by definition) is contained as a factor in 
j=l J 

several of these equations. 

For the first-order recombination at the surface, the following 

w 

relation between the atom mass flux and the wall number density n A 

may be employed 

% * - k wA < n A m A> < 33 > 

where the minus sign indicates that the flux is in the direction opposite to 
the outward normal and the expression for the rate constant k w ^ with 
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diffusion and slip is (Ref. 2). 


tt 



( 34 ) 


Here y a is the recombination coefficient. 

For a fully catalytic wall (y a = 1), the maximum value of the rate 
constant k wl . is limited by the surface temperature. The reaction rate con- 
stant for a fully catalytic wall with the gas phase in chemical equilibrium 
is often assumed to be infinity for the sake of simplicity. 

Employing Eq. (33), the net mass flux appearing in Eqs. (24) and 

(32) may be defined as 



“ k wA n A 



A = 0, N 


(35a) 


% = k wA n A m A ; 


r A = 0 for M = 0 2 
l A = N for M = N 2 


(35b) 


For all other atoms and molecules. 


M. = 0 


(35c) 


Equation (24), with M. defined by Eqs. (35a) through (35c), gives 


ffAs shown in ref. 2, by neglecting slip but keeping diffusion, a slightly 
different form of Eq. (34) is obtained: 


2y A, / kT w 


wA 


2-Y n 


2nm, 
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the number density ratio n*/n^. However, to obtain n? from this ratio, 

s 

an expression for n. is required. 

The net mass flux of 0 and N atoms to the surface, M Ay , is also 
equal to the rate of consumption of these atoms at the wall from surface 
recombination: 


% * *A M Ay : A * °» N 


(36a) 


and the corresponding net mass flux of Qfe and N 2 molecules will be 


M = . v • / A =* 0 for M =* O 2 

My y A Ay ’ l A = N for M * N 2 


(36b) 


For all other species, the net mass flux to the surface may be assumed to be 
zero. Thus, 


M iy*° 


(36c) 


Substituting values of the net and incident fluxes M.. , M. from Appendix 

y y 

B in Eq. (36), the following expression is obtained (after we neglect 
thermal diffusion): 


NS c 3C, NS 

l °ij [ c j l ( 

]=1 J 3y J q=l 


NS = 3C_ 

m q 


J s 
j*i 


m q 3 y 


•)] 






(37) 


IS 



which may be used to obtain 


s § 


Here 


is the source term defined as: 


* 


s _ y A 
A 


i yi 

VT 2S s 



c 



(38a) 


For Q 2 and tfe molecules 



(A = 0 for M = O 2 
1 A = N for M = N 2 


(38b) 


For all other species 

^ = 0 


(38c) 


SIMPLIFICATIONS FOR A MULTICOMPONENT MIXTURE 
Equations (24) (29), (30), (31), and (32) for multicomponent gas flows 
can be simplified if one makes the following assumptions: 


(i) All the diffusion coefficients, D. for a multicomponent gas have 

■ J 

^If no assumption is made about Djj, Eq. (37) would give an expression for 
(3 C.j/3 y) s for all the species: J 


3 C . m . 

(rrl, -(^ 


'ay 


m. 



NS S 

2 9 .], 

# m q ^ ' 


- 4 / 


m-i 


NS 

A (D U' C J^ 

3*1 


The source term in this expression, however, may be simplified to yield 

an expression for C| (or n|) for the recombining atoms only : 


r s m A n A 
u A - 


2-Y m. 

- 2 V¥ (_i) (_i) 


m. 


r”m7~ NS . 3 C . NS 
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+ 

the same value so that D. . - Di 2 , Di 2 is the same as the 

binary diffusion coefficient D l2 . 

(ii) The normal momentum flux to the pressure ratio, P^y/P-j is the 
same for all species and equal to that of the mixture. This also 
implies that the normal shear stress t . for species i is the 

J J * 

same as that for the mixture (t ). 

J J 

(iii) The rotational and vibrational states are fully excited so that 
the internal energy e^ for the air molecule may be taken as 
equal to 2kT/m.p 


*This is a somewhat stronger assumption. Because D-j -s are concentration de- 
pendent, whereas z? .s are virtually independent of composition. The multi- 
component diffusion J coefficient is related to the binary diffusion 

coefficient D ^ through the following relations (see Ref. 7): 


D ij a K ij “ (V M j) 


where quantities K-y are coefficients in a matrix which is the inverse of 
the matrix with the following coefficients: 


K 




+ Mj 


£1 
q*i 



(1*J) 


(1-J) 


One can see, therefore, that by employing D 12 (which is same as ^ 12 ) 
for all the species in a multicomponent gas mixture, considerable saving is 
obtained in computational effort and time without losing the general flavor 
of multicomponent diffusion. This is particularly true if the dissociated 
air consists predominantly of nitrogen molecules and oxygen atoms. The 
shuttle entry conditions fall into this category. Simple and multicompo- 
nent diffusion gave same results in "An Experimental and Analytical Study 
of Slip and Catalytic Boundary Conditions Applied to Spheres in Low 
Reynolds Number Arc Jet Flows," by Carl D. Scott, proceedings of the 9th 
International Symposium on Rarefied Gas Dynamics, Gottingen, July 15-20, 
1974, pp. 0.14-1 to 11. 



These assumptions retain the major effects of multicomponent fluxes on 
various slip quantities and provide considerable saving in computational 
effort required for the analysis of a flow problem (ref. 2 and footnote on 
page 20). With these simplifications, Eqs. (24), (29), (30), (31) and (32) 
yield: 

Concentration slip 
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Temperature slip 
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or, if the internal energy is frozen during reflection from the surface (see 
Eq. (28) of Ref. 5, for example). 
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Equation (37) may also be simplified to yield an explicit expression for no- 


where 


— ui), 


ac. 

( 1) + 

P. 3y S rn. D - 

± [ 1 + 1 — ii-.] 
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i p ■ #m q, 
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(45a) 
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s . 
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f A - 0 for M = 02 
1 A = N for M = N 2 


(45b) 


3 0 for all other species (45c) 

It is suggested here that the concentration for the major species (for 

example, nitrogen) be obtained by requiring the sum of concentrations of all 

\ 

the species to equal unity. It may be mentioned here that the mass of the 
ith species, m., is related to the molecular (or atomic) weight, W.., 
through the relation 

m. 


W 


i 


k. 

R 


(46) 


where k is the Boltzmann constant and R is the universal gas con- 
stant. 


CONCENTRATION SLIP BOUNDARY CONDITION FOR A FULLY CATALYTIC 
AND A NONCATALYTIC SURFACE 

s 

Equation (44) gives slip values of the concentration n^ for a finite 
catalytic surface. For a fully catalytic (y^ = 1) surface one generally 



assunes complete recombination of atoms at the surface. There is a slight 
discrepancy in this assumption because the maximun recombination rate is 
limited by the surface temperature as discussed earlier [See the discussion 
following Eq. (34)]. Thus, for a fully catalytic surface, equation (44) or 
(41) should be employed with y a = 1. 

For a noncatalytic surface (y a = 0), Eq. (45) gives = o for all 
the species. For this case, Eq. (44) becomes 
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which may also be written as 


(47a) 
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(47b) 


Summing the above equation over all the species gives 


NS = SC 

l (— 

q*l m q 3y 


(48) 


for a noncatalytic surface. Therefore, Eq. (47a) is not an appropriate 
boundary condition. However, employing Eq. (48) with Eq. (47b) yields: 



which may be used as the boundary condition for a noncatalytic surface with 
a multicomponent gas mixture. 

SIMPLIFICATIONS FOR A BINARY MIXTURE 


At lower altitudes with the flow in slight nonequil ibrium (i.e. when 
the Reynolds number parameter e 2 << 1), the derivatives of various quantities 
with respect to x and z and some other higher order terms like 3Vgy/3y 
may be dropped through an order of magnitude analysis (see Ref. 9, for exam- 
ple). Further, the dissociated air may be considered as a binary mixture 
(i.e. consisting of atoms and molecules only, see Ref. 8) at these altitu- 
des. If an assumption is also made that the internal energy of the mole- 
cules remains frozen during reflection from the surface, Eqs. (39) through 
(45) can be simplified further to the forms given here. 
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Velocity Sl ip 
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(52a) 


(53a) 


(54a) 



or, with frozen internal energy during surface reflection. 
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In Eqs. (54a) and (54b), n A is obtained from n A by using Eq. (50). 
Expression for n A 


n s , p s rr 1 . 1 1 

A 




3C, 


m 


A m M, A 


% % 3 y 


3C, 


-Y.C./JiL (^_) + _L3)] 


2irm A pp, 


'AM 


2m A 3y 


3C, 


/[(-*--!) (-*),] 
™A "m 3 * 


(55a) 


which may also be written as 



(55b) 
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A m A „ 
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For a fully catalytic surface (y^ = 1), Eq. (55b) gives appropriate 

t 

value for the concentration slip. Hie relevant boundary condition for a 
noncatalytic surface (y a * 0) will be 

3C « 

(— ■), - 0 (56a) 

3y 

as can readily be seen from Eq. (55b). 

For a noncatalytic surface (y a = 0), pressure and temperature slip Eqs. 
(51) and (54b) are further simplified, with the help of Eq. (56a), to: 
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Temperature slip (with frozen internal energy during reflection from the 
suface) 
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(58a) 


^Sometimes in the literature C A is prescribed as zero for a fully cata- 
lytic surface. Strictly speaking, this will be true only when the 
Reynolds number parameter e 2 (Ref. 8) is approximately zero (close to the 
chemical equilibrium condition at low altitudes). This can easily be seen 
by nondimensional i sing equation (55b) in a way similar to equation (28d) of 
reference 8. Thus, the recombination rate coefficient Y/\ and density (as 
measured through e 2 ), both control the recombination rate and not y A alone. 



Equations (52a), (53a), (57a), and (58a) ( with the noncatalytic surface 
assumption in Eqs. (57a) and (58a)) can be rewritten as: 
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Pressure slip 
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The concentration slip condition consistent with Eqs. (52b), (53b), 
(57b), and (58b) is 


(— ) . o (56b) 

ay 5 

The range of values of Ai , Bi , and Cl is as given below for a mixture of 
oxygen atoms and molecules and y 3 1.4 (ref. 8): 


1.0039 < Ai < 1.0186 

0.9507 < B l < (0.950 7)x(JI ^5. +• JT ^L) /( — + — ) 

^4 m A N 2 m M J l m, m, J 

0.9056 < Ci < 0.9507 


The minimum values for Ai and Ci occur at n A ■ 0.5 p /m^, whereas the 
maximum value for Bi occurs at this value of n^. 

The expression (52b), (57b), and (58b) reduce to those obtained in Ref. 
ft 

8 if one assumes the values for \ , B x , and C x to be unity and a non 

^See Appendix C for the dimensional form of the slip boundary conditions 
given in Ref. 8. 
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catalytic surface boundary condition as given by Eq. (56b). The concents 
tion slip boundary condition provided in Ref. 8 is: 


2-y i 


cl = c„ + r ill 

A Ae ^ J 
2y a 


2irm, 


kL 


AM 
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(59) 


which is the same as Eq. (55b) obtained here except for the additive term 
s ttt 

C^. It is obvious that the temperature slip boundary condition of Eq. 
(2.8b) obtained in Ref. 8 is valid, strictly speaking, for a noncatalytic 
surface only and is not consistent with the concentration slip boundary 
condition of Eq. (2.8c) (reproduced here as Eq. (59)) obtained for a finite 
catalytic wall. As a matter of fact, the velocity and temperature slip 
boundary conditions of Ref. 8 are similar to those of Ref. 9,^ where these 
slip conditions are provided for a perfect gas (or single species mixture). 

The inconsistencies in the boundary conditions used in Refs. 10 and 11 
are similar to those of ‘Ref. 8, namely, the pressure and temperature slip 


^It appears that Eq. (2.8d) of Ref. 8 for the concentration slip has been 
formulated for small deviations from the chemical equilibrium condition 
(i.e. the Reynolds number parameter e 2 « 1). Thus, when flow goes to 
chemical equilibrium (with e 2 * 0), one obtains from Eq. (2.8d) the 
equilibrium value for the concentration i.e. = C^, which would be 
zero for the oxygen atoms for surface temperatures of 2000 # K or less. 

^There appears to be some error with the form of Eqs. (2.7c) and (2.7d) 
given in Ref. 9 if one employs the definition of dimensionless heat-trans- 
fer rate * q ' In these equations from Eq. (2.4b) of the same reference. 

(The Reynolds number parameter, e 2 , is missing in Eqs. (2.2b) and (2.4b) 
for dimensionless shear stress and heat-transfer, respectively, of the same 
reference. These have been corrected in Ref. 8.) The resulting pressure 
and temperature slip expressions will contain the dimensionless viscosity 
coefficient 'm ' without the Prandtl number, see Ref. 10 and 11, for exam- 
ple. This is in contradiction to Eq. (2.8c) of Ref. 8, which contains 
dimensionless ’ K* in place of dimensionless V. For a perfect gas, 
employing 'y ' in place of ' K ' in the pressure and temperature slip equa- 
tions will result in no error. However, for real gas flows, the pressure 
and temperature slip values will be in error by the factor of y/K. The 
present analysis indicates that Eq. (2.8c) of Ref. 8 is of the "correct" 
form for a single species mixture. 



values are in error by a factor of dimensionless viscosity to the thermal 
conductivity ratio 'v/K'. References 10 and 11 also employ the single 
species formulation for obtaining the velocity, temperature and pressure 
slip values and the concentration slip expression in these references is 
that obtained for a binary mixture (Ref. 8). The coefficients appearing in 

the slip equations of Ref. 11 may be obtained by multiplying ai , bi , and 
c x (coefficients of Ref. 8) given here by 16/5jr. The factor 16/5 is some- 
times replaced by ir . 


SIMPLIFICATIONS FOR A SINGLE SPECIES MIXTURE 


Eqs. (39) through (42) and Eq. (43b) may be simplified for a single 
species mixture to the expressions obtained in Ref. 6. For small jump (or 
slip) conditions these simplified equations may be written as (with i A = 0 
and M.y * 0 for a single species mixture): 
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where we have used the following relations between the coefficients of vis- 
cosity and thermal conductivity and the mean free path for perfectly elastic 
spherical molecules (ref. 12). 
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The expression for u given earlier, following Eq. (25), for a single spe- 
cies mixture is modified to: 

u - i. nkTb 0 


If, now, an assumption is made that the slip values are small , equations for 
the density, pressure, and temperature slip are simplified to: 


* 4 * 

'Since the temperature slip expression contained in Ref. 6 is for a gas 
consisting of perfectly elastic spherical molecules, the gas does not 
possess any internal degree of freedom. Therefore, we have obtained the 
temperature slip equation from Eq. (43b). 
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Density slip 
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(64c) 


In obtaining Eq. (64c) from (64b), we have used Py/p S * 1 in the second- 
order terms (i.e. 2nd and 3rd terms on right side of Eq. (64b)). Equations 
(61b), (62), (63), and (64C) are the slip equations given in refs. 6 and 
12. 


NO-SLIP 80UNDARY CONDITION FOR THE SPECIES CONCENTRATION 
Multicomponent Mixture . The no-slip boundary condition for the species 
concentration without any assumptions may be obtained from Eqs. (44) and 



(45). In the absence of slip, the Knudsen layer thickness shrinks to zero 
and the values at the top of the Knudsen layer become the wall values (see 


Fig. 1): 


with 


ac. , 1|». 

c— 4 W ♦ i- (— ) w 

w p w r 3y m i d 12 , 

n? - -1 [ 1 + 2 — ] 


■1 


m 


i 


NS _ 3C_ 

l (*- — 3.) w 

q =l ^ 3 y 


( 66 ) 


w 


2y / 


1 m A n A 


(2-y.) * 2 n 



A = 0, N 


(67a) 


^w = , ^w . rA = 0 for M ■ 0 2 . 

V A * l A = N for M = N 2 * 


(67b) 


,w 


ij>.j = 0 for all other species 
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w 

where we have neglected the higher order shear (i.e. P^/p - 1). 

It is suggested here, again, that the concentration for the major 
specie (for example, nitrogen) be obtained by requiring the sum of concen- 
trations of all the species to equal unity. 

For the recombining 0 and N atoms, Eq. (66) may further be simpli 
fied to: 
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The recombination rate constant k wA in Eq. (68) has been defined as (Ref. 
2 ) 
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fJ * incident flux 

f* » specularly reflected flux 

pw * diffusely reflected flux 



TCy) 


Figure 1. The Knudsen layer showing general fluxes and coordinate axes. 

The temperature as a function of normal distance is schematically 
overlayed. 
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without slip and higher order shear. 

For a none atalytic wall (with y a = 0), Eq. (49) gives 
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for all the species of a multicomponent mixture. 

Binary Mixture . For a two-species mixture of atoms and molecules, Eq 
(55b) gives 
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for a surface with finite catalycity and 
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for a noncatalytic surface. Eqs. (71) and (72) are employed in Ref. 4. 

Appendices D and E give the slip and no-slip boundary conditions (pre- 
sented in the text earlier), in the dimensionless form for the body-fitted 
and spherical polar coordinates, respectively, for a planar flow. The 
various integrals employed for evaluation the net, incident, and specularly 
reflected fluxes defined through Eqs. (5) to (7), respectively, are given in 
Appendix F. 
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DISCUSSION AND CONCLUSIONS 


The present analysis provides a consistent formulation for the slip 
equations for a multicomponent, binary, and single species mixture reacting 
catalytically on the surface. The differences between the slip expressions 
obtained by various researchers have been reconciled and implication of 
various assumptions (some of them inconsistent) contained in those express- 
ions is discussed. The slip equations have also been obtained in body-fit- 
ted and spherical polar coordinates in a form which can readily be employed 
in the flow-field calculations. 

Usually, the equation for pressure slip is not required as a boundary 
condition, but is needed to obtain the surface pressure. The temperature 
slip equation given here is for a constant surface temperature, which is 
provided as a boundary condition. For an adiabatic surface, however, the 
slip temperature, T s , may be obtained by equating the wall heat-transfer 
rate to zero 1 .e. 
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■) s = 0 


where the expression for mass flux is provided in Appendix B and the 
higher order terms have been dropped. The temperature slip equation will 
now be required to obtain the wall temperature, T . 

An expression has also been obtained for the finite-rate species-con- 
centration boundary condition for a multicomponent gas mixture without sur- 
face slip. This boundary condition in the literature (Ref. 4) has generally 
been specified by assuming the dissociated air as a binary mixture. 



However, the binary mixture formulation does not provide boundary conditions 
for the recombined molecules (on surface) such as Ofc and N 2 in a multicompo- 
nent mixture; it gives boundary conditions only for the recombining atoms. 

To assess the importance of various terms in the general boundary con- 
ditions suggested here, a study is being undertaken in which the flow prop- 
erties and boundary conditions would be varied systematically. The boundary 
equations form a simultaneous set, which is being coupled with a flow field 
calculation procedure in the stagnation region. These equations would 
finally be coupled with the viscous shock-layer code developed by Moss (Ref. 
13) for the detailed flow-field calculations. The boundary equations 
obtained in the present work should provide a more realistic set of boundary 
conditions for a multi -component mixture for low-Reynolds-n umber slip flows 
as well as no-slip flows. 

In conclusion, the boundary slip expressions obtained here are closed 
form solutions of the mass, momentum, and energy flux equations using the 
Chapman-Enskog velocity distribution function. This function represents a 
solution of the Boltzmann equation in the Navier-Stokes approximation. 
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APPENDIX A 


CONSTANTS APPEARING IN THE DISTRIBUTION FUNCTION 
AND DIFFUSION VELOCITY IN TERMS OF TRANSPORT PROPERTIES 

t * i 

The Sonine expansion coefficients a^, a.^, b.g, and c^q appearing 
in the general velocity distribution function are found by a variational 
technique in which they are solutions to sets of simultaneous equations. 
References 2 and 7 provide the coefficients in terms of solutions to these 
set of equations. These solutions are expressed in terms of the collision 
integrals . Reference 7 also provides the transport properties in 

terms of the Sonine expansion coefficients. Thus, in place of evaluating 
these coefficients in terms of the collision integrals, they may be 
expressed in terms of the transport properties. The various relations 
are: 


( i) for a. 0 



(Al) 


(ii) for a n 


K(0 * 



(A2) 


(iii) for b iQ 

. NS 

u(«) - - kT I n. b, 0 (5) (A3) 

2 i=l 1 1U 

(iv) for c^ 

In the distribution function used here the Kernel 



(A4) 


D ij(5) 



( 5 ) 


where dI (5), K(5),u(5), and D i - (5 ) are the multi -component thermal 

diffusion coefficient, thermal conductivity, coefficient of viscosity, and 
the multicomponent diffusion coefficient, respectively. The argument 5 is 
the number of terms used in the Sonine expansion. Except for dT, letting 
5=1 gives quite good results for K, u and D... When 5=1, however, 

' J 
T 

coefficients vanish. Hence, in order to get the coefficient of thermal 
diffusion, it is necessary to take at least two terms in the Sonine expan- 
sion (i.e. 5 = 2). If the argunent does not appear with a coefficient (ex- 
cept for a.j 0 ), it is considered to have one term in the Sonine expansion. 

The diffusion velocity for ith species in terms of the transport co- 
efficients is obtained as (see ref. 7) 


- MC 

*1 - (-$ i i - — — 

n^p j=l J ,J * n i m i 1 3x k 


(A5) 


where the diffusion vector d k for the jth species is defined after ne- 
glecting the external forces as (see Eq. (16) in the text): 


„ n . n . n .m . 

d J = JL_ ( 4 ) + ( 4 - - -i-i) — (tnp) 
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Eq. (A6) may be simplified to obtain 


, 2 i 3C, NS , 3C n 

J. C( i (J. -i)] 

n 2 "j S 3x k 


J I ' 1 ™q 3x k 


(A7) 


SIMPLIFIED EXPRESSIONS FOR COEFFICIENTS a n AND b iQ 
Since a.j and b.g can not be obtained directly from equations (A2) 
and (A3), we introduce the following assumptions for the mixture thermal 
conductivity and viscosity; respectively. 


NS n. 

I 4 K n * 

i=l n 1 


(A8) 


NS n. 
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i 3 l n 
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Equations (A8) and (A9) are approximate forms for the more exact formula of 
Wilke (see ref. 14). These equations imply that 


NS n. 

I 4 

j=l n 


ij 


(A10) 


in the Wilke's relation. Thi s is approximately true for air. 

Thus, with the help of Eqs. (A2), (A3), (A8), and (A9) we may obtain 
the following approximate expressions for a^ and b^, respectively. 


n i a il " ' 


. n, k. isr 
1. 1 1 / 1 
5 n k 
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n i b iO 
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kT n 


( A1 2) 


For accurate evaluations of a^ and b. Q , cumbersome expressions of the 

type of Eq. (7.4-49) of ref. 7 are required to be solved. 

/ 



APPENDIX B 


EVALUATED EXPRESSIONS FOR THE FLUXES OF MASS, 

MOMENTUM, AND ENERGY FROM INTEGRALS OVER 
THE DISTRIBUTION FUNCTION 

The fluxes of mass, momentum, and energy (i.e. of a property 
are given in terms of integrals over the velocity distribution function 
f(V) in the "Analysis" section of main text. Here we give evaluated forms 
for these fluxes obtained from integration over the distribution functions. 
The various integrals needed in these evaluations are provided in Appendix 
F. 

(i) Expressions for Mass Flux of Species 
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Specularly Reflected 



Diffusely Reflected 


1 2yT 


2 k T 1/2 


m. 


i 


(B4) 


It may be mentioned that as obtained here is similar to the mass flux 

(j^) expression of Ref. |£ 


(ii) Correlation with Blottner's Expression (Ref. 15) for Mass Flux 

If the thermal diffusion term is neglected in Eq. (Bl) and 
substitutions are made for c^ and d^ from Eqs. (A4) and (A7), Eq. (Bl) 
may be simplified to obtain* 
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where is the molecular (or atomic) weight of the ith species and is 
related to the mass of the ith species, m^ , through the relation 



(B6) 


Here k is the Boltzmann constant and R is the universal gas constant. 
Equation (B5) may also be written as 
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or, with the introduction of the multicomponent Lewis nunber, L.., defined 

I J 


*It may be noted here that 0. See reference 7 for details. 
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Eq. (B7a) may be rewritten as 
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where 


In general , 
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^See Eqs. (7), (8a), and (8b) of reference 15. 



If the binary Lewis numbers are assumed to be the same for all the 

species, then 


L.. * L.. = L- 4 (C. , M., M.) 
ij ij l* i* j' 


and Eq. (B8) for Ab^ becomes 
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and the mass flux Eq. (B7e) now becomes 
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Further, if the binary Lewis numbers, L.y are constant for all the 
species, the term 
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is zero in Eq. (BIO). Therefore, 
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Thus, Eq. (BIO) gives the mass flux due to concentration gradients only with 

•Jt 

and Le^ defined through Eqs. (B9) and (Bllb), respectively. 

★ 

To utilize Eqs. (B9) and (Bllb), are still to be evaluated 

through complex matrix inversions (ref. 15). If we now make a little 

★ 

stronger assumption such that Lj • = L i2 a Le ( same f° r the species). 



we obtain a much simpler expression for Ab^. 

and 

Eq . 
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Similarly the expression for Ab^. may be obtained from Eq. (B9) as 
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It may be noted that for the case when L = 0, the term Le. — 1 in Eq. 

ay 3y 

(BIO) vanishes and Eq. (B12b) is not required. Further, employing a 
constant value for the Lewis number, Le (= p Pr £ 12 /y), does not imply Le. 
to have a constant value. 

It should be pointed out here that for the case when L-. is the same 

'J 
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for all the species we have used the condition 


NS 
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j*i 
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3 ♦ 

(and not just Ab ik = 0 as employed in ref. 15) to obtain Eq. ( 81 lb) for 
Le.. . The present formulation appears to provide a more consistent 
expression for Le. . 

Finally, Eqs. (BIO), (B12b), and (B13) for a binary (two species) 
mixture give 

3C. 

j. * — L Le ( — ) ; (Fick's law for mass flux) (B14) 
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(iii) Expressions For Normal Momentum Flux 
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where are terms in the stress tensor for the ith species. 
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(iv) Expressions For Tangential Momentum Flux 
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(v) Expressions For the Energy Flux 
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APPENDIX C 


SLIP BOUNDARY CONDITIONS GIVEN BY DAVIS FOR A BINARY MIXTURE 
In this appendix, the dimensional form of the slip boundary conditions 
obtained by Davis (Refs. 8 and 9) is provided for the purpose of comparison 
with the slip expressions obtained in the main text for a binary mixture. 

The dimensionless form of the surface slip conditions provided in refs. 
8 and 9 are 

Concentration slip (ref. 8)* 



Velocity slip (ref. 8) 



(C2) 


Pressure slip (ref. 9) 


tt 



*Here the dimensionless quantities are denoted with a 'bar*. See Appendix 
D for details of nondimensional ization and the coordinate system. 

++ — . 

For the reasons explained in the main text, we have employed 'K' in place 
of IT' in this expression. 



Temperature slip (ref. 8) 
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With the help of the nondimen si on al i zing quantities given 
the following dimensional forms of Eqs. (Cl) through (C4) 
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in refs. 8 and 9, 
are obtained: 
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* 

where n is the coordinate normal to surface (see Fig. D1 in Appendix D) 

33 = 

and n and m are the number density and mass of a mixture molecule, 
respectively. The constants ai , bi , and ci are defined as: 
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APPENDIX D 


SLIP CONDITIONS IN THE BODY-ORIENTED COORDINATE SYSTEM 
In this Appendix we obtain the slip conditions for an axisymmetric body 

★ Hr 

in the body-oriented coordinate (s , n ) system. The coordinate configura- 
tion for such a system is given in Fig. Dl. The equations of this Appendix 
contain simplifications similar to those given for a multicomponent mixture, 
following Eq. (38c) of the main text. 



Figure Dl. Coordinate Configuration For A Body (Symmetric About Its Axis) 

By employing the metric coefficients (see ref. 16) of h x = 1 + n < and 
hz 3 1, equations (39) through (45) give the following slip equation in the 
body-oriented coordinate system. 
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If the internal energy (comprising of rotational and vibrational energies in 
this case) is considered frozen during reflection from the surface, the 
following form of the temperature slip equation is obtained (see Eq. (28) of 
Ref. 5): 
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In Eqs. (Dl), (D5), (D6) and (D8a) the ratio ~ is defined as 
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ana b Q is related to viscosity, u, through the relation 
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which has been obtained from (A3) by assuming b^ to be the same (bQ) for 
all the species. 

Equations (Dl) through (D8) may now be nondimensional ized by using the 
following relations (refs. 8 and 13): 
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Introducing the nondimens ional quantities as defined here, the following 
equations are obtained from Eqs. (D-l) through (D-8): 
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(012b) 
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Or, with frozen internal energy during reflection from surface. 
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(Prandtl lumber) 
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(Recombination rate 
constant) 


For a noncatalytic surface (y a ■ ky A 3 0), eqs. (D18a) through (D18c) give 
^ = 0 for all species. For this case, Eq. (D17) becomes 
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NS U 3 C 
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3n 
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1 q*l Wq 3n s 


Equation (D20), similar to Eq. (47b), gives 


(D20) 



(D21 ) 


which may be employed as the boundary condition for a noncatalytic surface 
with a multicomponent gas mixture (with the binary assumption for diffusion 



coefficients). For a fully catalytic surface, employing y a 3 1 and using 
Eqs. (D18a) through (D18c) in Eq. (D17) would yield the appropriate 
concentration slip. 

If the multicomponent diffusion coefficients are retained without the 
binary assumption the underlined terms in Eq. (013) are to be replaced by 

NS l,, 3C, NS a 9C. 

- I -ii [ _JL - c, l L -1] ■ 

j a l Le 37T J q=l W„ 3rf s 

j*i q 

and the underlined termsin Eq. (D14) are to be replaced by 
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where the multicomponent Lewis nunber, L. is defined as 


*1j 
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pC D. . 
P U 

K 


Simplifications For A Binary Mixture 

When all the species in a gas mixture can be considered as atoms and mole- 
cules only (see ref. 8), Eqs. (Dll) through (D18) may be further simplified. 
With the assumption of Py/p S " 1, Eqs. (Dll) through (D18) for a binary mix- 
ture yield: 



Concentration Slip 



Pressure Slip 


(D22) 



4 

5 VST 


( JLJ (i±) iL 

Y-l s 8 Pr 



W 

x(_f}3/2 

W A 
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where we have used the equation of state 



1 r T s T refj ± 

M 2 - T m T. fl s 


(024) 


and neglected the higher order shear terms as well as diffusion terms. 


Velocity slip 




k u 




l+n < 


(D25) 


where we have again employed the equation of state given earlier and neglec 
ted higher order shear, conduction, and diffusion terms. In obtaining 
Eq. (D25) we have also used the approximation: 
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X 


2 


1 


[C A S (W M -W A ) +„/'» 

■ « 

[ C A W^M ” +y ^p) 


To be more exact one may keep this factor in Eq. (D25). 


Temperature slip 



(026) 


where 


[ C A (W M - H A ) + w a ] 3/2 

3/2 3/2 3/2 

[C A (W M - 2 W A ) + 2 W A ] 

and ~ < xi < 1 for 0 < C A < 1 and diatomic molecules (W M * 2 W A ) . In 

obtaining Eq. ( 026) , we have employed the equation of state and = 2W A . 

With frozen internal energy during reflection from the surface, the 
temperature slip equation becomes 
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where 


X 2 


CC A ( - «A> * « A ] 3/2 
r , 3/2 3/2 3/2 

[ C A (W M - W A ) + W A | 


For diatomic molecules (W M * 2 W A ), X 2 ma y be taken as unity for 
0 < C A < 1. Equation (D27) employs = 2W A< 


s 

Equation for C A 



(D28) 


Equation (28), for a noncatalytic surface with y a 3 k wA 3 0, gives 
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3“n s 


(D29) 
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whereas using = 1 in Eq. (D29) and in the expression for k wA : 


wA 
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T 
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w 
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(D30) 


would yield the appropriate concentration slip for a fully catalytic sur- 
face. 

With the following {somewhat inconsistent) assunptions, Eqs. (023), 
(D25), (D27), and (028) may be simplified to those obtained in refs. 8 and 
11: 

(i) In pressure slip Eq. (023), Pr = 1 is employed along with approxi 
mations 



(ii) In velocity slip Eq. (D25) the following assumption is made for the 
mixture molecular weights 
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(iii) In temperature slip Eq. (D27) an assumption is made for the surface 
to be noncatalytic so that k w ^ * 0. In addition, Prandtl number 
(Pr) is taken unity and it is assumed that 


1 

2 



15 

8 


(iv) No assumption is required in the concentration slip Eq. (D28) . But 
it may not be appropriate to employ (Ref. 11) this equation for all 
the species of a multicomponent mixture. This equation is applic- 
able to concentration slip of recombining atoms only in a binary 
(two-species) mixture (Ref. 8) of atoms and molecules. 

No equation has been obtained in reference 8 or 11 to correspond to Eq. 

w 

(D22) to obtain wall values of the species concentration (C^) from th§ 
values at the edge of the Knudsen layer (C^). 


Slip Expressions For A Single Species Mixture 
For a single species mixture (y a = 0), the following slip expressions are 
obtained: 



(031) 
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Pressure Slip 
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Velocity Slip 
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Temperature Slip [For a gas consisting of molecules only (i.e. diatomic 
perfect gas)] 
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where we have used xi a — » or with frozen internal energy during 

2 

reflection from the surface. 
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S W 2 N 2 y-1 0 
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s s 
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Equations (D32) (D33) and (D35) are the ones employed in References 8 and 11 
with Pr ■ 1. . 


No-Slip Species Concentration Boundary Condition 
Multicomponent Mixture 

The no-slip boundary condition may be obtained from Eqs. (D17) and 
(D18) . in the absence of slip, the Knudsen layer thickness shrinks to 
almost zero, the values at the top of the Knudsen layer become the wall 
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values (See Fig. 1) : 
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where we have neglected the higher order shear (i.e. Py/p - 1). 
For the recombining 0 and N atoms, Eq. (D17) may also be 
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(D36) 


(037a) 

(D37b) 

(D37c) 

written as 
(D38) 


The recombination rate constant k wA in Eq. (D37a) and (D38) has been 
defined as (Ref. 2) 



(D39) 



( *Aj J_ f *~ T w 
2-Y a M, 1*7 m W A 


without slip and higher order shear. 

For a noncatalytic wall (with y a 3 k wA =0), Eq. (D21) gives 
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3C 


3TT 


i) 
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(D40) 


for all the species of a multicomponent mixture with the binary assumption 
about the diffusion coefficients. 

Binary Mixture 

For a two-species mixture of atoms and molecues, Eq. (D28) gives 



(041a) 


which may also be obtained from Eq. (038) for a surface with finite 
catalycity. k wA in Eq. (41a) is again obtained from Eq. (039). For a 
noncatalytic surface with IT . * 0, Eq. (D41) gives 


3C a 

i-A v ■ 0 

3n w 


(D42) 
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Eqs. (D41) and (042) are obtained in ref. 4. However, ref. 4 employed Eq. 
(D30) in place of Eq. (039) for obtaining the reaction rate constant, k wA 
As obtained in ref. 2, Eq. (D30) Is more appropriate when slip and higher 
order shear are included in the flowfield analysis (e.g. under rarefied or 
low density conditions). 

If k wA is substituted from Eq. (D39) in Eq. (041a), we may also 
obtain 



where we have used the equation of state 



_1 

OP 





(D41b) 


(D43) 


Eq. (D41b) compares with the corresponding equation of ref. 8, if one 
keeps in mind that the diffusion coefficient F Q (of Ref. 8) in the absence 
of slip is related to the present variables by the relation 



- , Lev 

W ' TT"7 w 
w pp w 


It may be seen from Eqs. (41a) or (41b) that the gradient (3C A /3n) w is 
governed by the ratio "k^/e 2 or y a /e 2 . Therefore, for surfaces which are 
almost noncatalytic, this ratio would be of the order of one for large 
values of the Reynolds number parameter (1/e 2 ). This would imply that a 



surface, regardless of its catalycity, would produce a larger effect 
on the concentration gradient for high density rather than low density 
conditions. 



APPENDIX E 

SLIP CONDITIONS IN THE SPHERICAL COORDINATES 
For analyzing the flow in the stagnation region of a body, it is conven 
ient to use the spherical polar coordinates (r,$) for the two-dimensional 
flow. The following relations exist between the spherical polar and the 

"At ★ 

body-oriented (s , n ) coordinates over the spherical portion of a body 





or, in the non-dimensional form. 


r sin $ = R + n cos <P 


r = r N + n = 1 + n ; k * 1 y 


<J> = S ; 3r = 3n ; 3<|> = 3 s 


(E2) 


where we have used nose radius, r N , to nondimensionalize all the distances 
(see Appendix D). 

Using the relations given by (E2) in Eqs. (Oil) through (D19), we can 
obtain the slip equations in the spherical polar coordinates for a multicom- 
ponent mixture. Once again, these equations employ the simplifications 
given for a multicomponent mixture, following Eq. (38c) of the main text. 
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or, with frozen internal energy during reflection from surface. 
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In obtaining these various equations, the following form of the equation of 
state has been employed: 


1 /■ T s ^"ref ^ r W< » ^ 


(E12) 


Similar to Appendix 0, Eqs. (E3) through (Ell) may further be simp! i> 
fied for a binary mixture and for a single species mixture. 
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APPENDIX F 
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INTEGRALS REQUIRED FOR OBTAINING VARIOUS FLUXES 
In order to obtain the normal fluxes of mass, momentum, and energy, it 
is necessary to evaluate various integrals over the velocity space of the 
distribution function. These integrals involve terms that are various 
velocity moments of the distribution function. The integrals are provided 
for the net, incident, and specularly reflected fluxes. These integrals are 
basically the same as those provided in Reference 2 except for some 
corrections and additions. 

FI. Integrals Over the Entire Velocity Space (For Net Quantities) 
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F2. Integrals Over the Lower Half Velocity Space (For Incident Quantities) 
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F3. Integrals Over the Upper Half Velocity Space (For Specularly 
Reflected Quantities) 
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(-W ) appearing in the integrals emphasizes that the sign of the thermal 
velocity component normal to the surface, V y (and consequently that of W y ) 
changes in the distribution function f upon reflection from the surface. 




89 



1. Report No. 

NASA TM-85820 


2. Government Accession No. 


4. Title and Subtitle 

Surface-Slip Equations for Multicomponent Nonequilibrium 
Air Flow 


7. Author(s) 

R. N. Gupta*, C. D. Scott**, and J. N. Mosst 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 


3. Recipient’s Catalog No. 


S. Report Date 
March 1985 


6. Performing Organization Code 

506-51-13-03 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


National Aeronautics and Space Administration 
Washington, DC' 20546 


15. Supplementary Notes 

Part of this work was presented as AIAA Paper 84-1732 at the AIAA 19th Thermophysics 
Conference, June 25-28, 1984, Snowmass, Colorado. 

*01d Dominion University, Norfolk, VA; **NASA Johnson Space Center, Houston, TX; 

JL.XT SO* T _ . _ 1 T-» 1 /-I . YT tT» 


I llSJ-kVlU. U V. J VC* 1. Vll V * — ' * * V V *• J * * V**** U V ** ) 

16. Abstract 

Equations are presented for the surface-slip (or jump) values of species concentra- 
tion, pressure, velocity, and temperature in the low-Reynolds number, high-altitude 
flight regime of a space vehicle. These are obtained from closed form solutions of 
the mass, momentum, and energy flux equations using the Chapman-Enskog velocity 
distribution function. This function represents a solution of the Boltzmann equa- 
tion in the Navier-Stokes approximation. The analysis, obtained for nonequilibrium 
multicomponent air flow, includes the finite-rate surface catalytic recombination 
and changes in the internal energy during reflection from the surface. Expressions 
for the various slip quantities have been obtained in a form which can readily be 
employed in flowfield computations. A consistent set of equations is provided for 
multicomponent, binary, and single species mixtures. Expression is also provided 
for the finite-rate, species-concentration boundary condition for a multicomponent 
mixture in absence of slip. 
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